Termination of Reentry in an Inhomogeneous Ring of Model Cardiac Cells 
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Reentrant waves propagating in a ring or annulus of excitable media are a model of the basic mechanism underlying a 
major class of irregular cardiac rhythms known as anatomical reentry. Such reentrant waves are terminated by rapid electrical 
stimulation (pacing) from an implantable device. Because the mechanisms of such termination are poorly understood, we 
study pacing of anatomical reentry in a one-dimensional ring of model cardiac cells. For realistic off-circuit pacing, our model- 
independent results suggest that circuit inhomogeneities, and the electrophysiological dynamical changes they introduce, may 
be essential for terminating reentry in some cases. 

PACS numbers: 87.19. Hh, 05.45.Gg, 05.45.-a, 87.10. +e 



Reentrant tachycardias, abnormally rapid excitations 
of the heart that result from an impulse that rotates 
around an inexcitable obstacle ( "anatomical reentry" ) 
[1] or within a region of cardiac tissue that is excitable 
in its entirety ("functional reentry") [2,3], may be fatal 
when they arise in the heart's ventricles. Trains of local 
electrical stimuli are widely used to restore normal wave 
propagation in the heart during tachycardia. Such "an- 
titachycardia pacing" is not always successful and may 
inadvertently cause tolerated tachycardias to degenerate 
to more rapid and threatening spatiotemporally irregular 
cardiac activity such as ventricular fibrillation [4]. The 
underlying mechanisms governing the success or failure 
of antitachycardia pacing algorithms are not yet clear. 
Understanding these mechanisms is essential, as a better 
knowledge of the processes involved in the suppression of 
ventricular tachycardia (VT) through such pacing might 
aid in the design of more effective therapies. 

Several factors influence the ability of rapid pacing to 
interact with VT. The most prominent are [5]: (a) VT 
rate (for anatomical reentry this is determined by the 
length of the VT circuit and impulse conduction velocity 
around the obstacle), (b) the refractory period (i.e., the 
duration of time following excitation during which car- 
diac tissue cannot be re-excited) at the pacing site and 
in the VT circuit, (c) the conduction time from the pac- 
ing site to the VT circuit, and (d) the duration of the 
excitable gap (the region of excitable tissue in the VT 
circuit between the front and refractory tail of the reen- 
trant wave [6]). A single stimulus is rarely sufficient to 
satisfy the large number of conditions for successfully ter- 
minating reentry. Therefore, in practice, multiple stimuli 
are often used - where the earlier stimuli are believed to 
"peel back" refractoriness to allow the subsequent stimuli 
to enter the circuit earlier than was possible with only a 
single stimulus [5]. 

The dynamics of pacing termination of one- 
dimensional reentry has been investigated in a number 



of studies [2,7-10], but most of these were concerned ex- 
clusively with homogeneous ring of cardiac cells. The 
termination of reentry in such a geometry (which is effec- 
tively that of the reentry circuit immediately surround- 
ing an anatomical obstacle) occurs in the following man- 
ner. Each stimulus splits into two branches that travel 
in opposite directions around the reentry circuit. The 
retrograde branch (proceeding opposite to the direction 
of the existing reentrant wave) ultimately collides with 
the reentrant wave, causing mutual annihilation. The an- 
terograde branch (proceeding in the same direction as the 
reentrant wave) can, depending on the timing of the stim- 
ulation, lead to resetting, where the anterograde wave be- 
comes a new reentrant wave, or termination of reentry, 
where the anterograde wave is blocked by the refractory 
tail of the original reentrant wave. From continuity ar- 
guments, it can be shown that there exists a range of 
stimuli phases and amplitudes that leads to successful 
reentry termination [7]. Unfortunately, the argument is 
essentially applicable only to a ID ring - the process is 
crucially dependent on the fact that the pacing site is on 
the reentry circuit itself. However, in reality, the location 
of the reentry circuit is typically not known when an elec- 
trical pacing device is implanted, and, it is unlikely that 
the pacing site will be so fortuitously located. 

In this study, we examine the dynamics of pacing from 
a site located some distance away from the reentry cir- 
cuit. Such off-circuit pacing introduces the realistic prop- 
agation of the stimulus from the pacing site to the reentry 
circuit. Because reentrant waves propagate outwardly 
from, in addition to around, the circuit, the stimulus 
will be blocked before it reaches the circuit under most 
circumstances. Multiple stimuli are necessary to "peel 
back" refractory tissue incrementally until one success- 
fully arrives at the reentry circuit [5]. Further, once a 
stimulus does reach the circuit, its anterograde branch 
must be blocked by the refractory tail of the reentrant 
wave, or resetting will occur and termination will fail. 
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However, as outlined below, this is extremely unlikely to 
happen in a homogeneous medium. 

Let us consider a reentrant circuit as a ID ring of 
length L with separate entrance and exit sidebranchcs 
(Fig. 1). This arrangement is an abstraction of the spa- 
tial geometry involved in anatomical reentry [5]. Further, 
let the pacing site be located on the entrance sidebranch 
at a distance z from the circuit. We use the entrance 
sidebranch as the point of spatial origin (x = 0) to define 
the location of the wave on the ring. The conduction 
velocity and refractory period at a location a distance 
x away (in the clockwise direction) from the origin arc 
denoted by c(x) and r(x), respectively. 

For a homogeneous medium, c{x) = c, r{x) = r (c, r 
are constants). Therefore, the length of the region in 
the ring which is refractory at a given instant is I = cr. 
For sustained reentry to occur, an excitable gap must 
exist (i.e., L > c r). We assume that restitution effects 
(i.e., the variation of the action potential duration as a 
function of the recovery time) can be neglected. Further, 
the circuit length L is considered to be large enough so 
that the reentrant activity is simply periodic. For conve- 
nience, associate t — with the time when the reentrant 
wavefront is at x = (i.e., the entrance sidebranch) [Fig. 
1(a)]. Let us assume that a stimulus is applied at t = 0. 
This stimulus will collide with the branch of the reentrant 
wave propagating out through the entrance sidebranch 
at t = z/2c [Fig. 1(b)]. The pacing site will recover at 
t = r and if another stimulus is applied immediately it 
will reach the reentry circuit at t = r + (z/c) [Fig. 1(c)]. 
By this time the refractory tail of the reentrant wave will 
be at a distance x = z away from the entrance sidebranch 
and the anterograde branch of the stimulus will not be 
blocked. Thus, when z > 0, it is impossible for the stim- 
ulus to catch up to the refractory tail in a homogeneous 
medium. This results in resetting of the reentrant wave 
rather than its termination. 

Note that, if the first stimulus is given at a time 
t < — z/c, it reaches the reentry circuit before the arrival 
of the reentrant wave. As a result, the retrograde branch 
collides with the oncoming reentrant wave, while the an- 
terograde branch proceeds to become the reset reentrant 
wave. Even in the very special circumstance that the 
reentrant wave reaches the entrance sidebranch exactly at 
the same instant that the stimulus reaches the circuit, the 
two colliding waves allow propagation to continue along 
the reentrant circuit through local depolarization at the 
collision site. As a result, pacing termination seems all 
but impossible in a homogeneous reentry circuit. 

The situation changes, however, if an inhomogeneity 
(e.g., a zone of slow conduction) exists in the circuit [Fig. 
1(d)]. In this case, the above argument no longer holds 
because the inhomogeneity alters the electrophysiological 
dynamics (notably refractory period) of the excitation 
waves. As a result, stimuli may arrive at the circuit from 
the pacing site and encounter a region that is still refrac- 
tory. This leads to successful block of the anterograde 
branch of the stimulus, while the retrograde branch anni- 



hilates the reentrant wave (as in the homogeneous case), 
resulting in successful termination. 

The following simulation results also support the con- 
clusion that the existence of inhomogeneity in the reentry 
circuit is essential for pacing termination of VT. We as- 
sume that the cardiac impulse propagates in a continuous 
one-dimensional ring of tissue (ignoring the microscopic 
cell structure) with ring length L, representing a closed 
pathway of circus movement around an anatomical ob- 
stacle (e.g., scar tissue). The propagation is described by 
the partial differential equation: 

dV/dt = -I ion /C m + DV 2 V, (1) 

where V (mV) is the membrane potential, C m = 1 fiF 
cm~ 2 is the membrane capacitance, D (cm 2 s _1 ) is the 
diffusion constant and I ion ([iA cm~ 2 ) is the cellular 
transmembrane ionic current density. We used the Luo- 
Rudy I action potential model [11], in which Ii on — I^ a + 
hi + Ik + Iki + Ik p + h- iNa = G Na m 3 hj(V - E Na ) 
is the fast inward Na + current, I S i = G S idf(V — E si ) 
is the slow inward current, Ik = GkxXi(V — Ek) is 
the slow outward time-dependent K + current, Iki = 
GkiK1oo(V — Eki) is the time-independent K + cur- 
rent, Ik p = 0.0183K p (V — Ek p ) is the plateau K + cur- 
rent, and I b = 0.03921(1/ + 59.87) is the total back- 
ground current, m, h, j, d, /, x and the gat- 
ing variables satisfying differential equations of the type: 
dy/dt = (yoo — u)/t v , where yoo and t v are dimensionless 
quantities which are functions solely of V . The exter- 
nal K+ concentration is set to be [K] = 5.4mM, while 
the intracellular Ca 2+ concentration obeys d{C&]i/dt = 
-1Q- A I SI +0.07(10" 4 - [Ca]i). The details of the expres- 
sions and the values used for the constants can be found 
in Ref. [11]. We solve the model by using a forward-Euler 
integration scheme. We discretize the system on a grid 
of points in space with spacing Sx = 0.01 cm (which is 
comparable to the length of an actual cardiac cell) and 
use the standard three-point difference stencil for the 1- 
D Laplacian. The spatial grid consists of a linear lattice 
with L points; in this study we have used L = 2500. 
The integration time step used in our simulations is St 
= 0.005 ms. The initial condition is a stimulated wave 
at some point in the medium with transient conduction 
block on one side to permit wavefront propagation in a 
single direction only. Extrastimuli are introduced from a 
pacing site. 

To understand the process of inhomogeneity-mediated 
termination we introduced a zone of slow conduction in 
the ring. Slow conduction has been observed in cardiac 
tissue under experimental ischemic conditions (lack of 
oxygen to the tissue) [12] and the conduction velocity in 
affected regions is often as low as 10% of the normal prop- 
agation speed in the ventricle [13]. This phenomenon can 
be reasonably attributed to a high degree of cellular un- 
coupling [14], as demonstrated by model simulations [15]. 
In our model, slow conduction was implemented by vary- 
ing the diffusion constant D' from the value of D used 
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for the remainder of the ring. The length and diffusion 
constant(D') of the zone were varied to examine their ef- 
fect on the propagation of the anterograde branch of the 
stimulus. We found that varying the length of the zone 
of slow conduction (specifically, between 1.5 mm and 25 
mm) had no qualitative effect on the results. For the sim- 
ulation results reported below we used D = 0.556 cm 2 /s 
and D' — 0.061 cm 2 /s, corresponding to conduction ve- 
locities c ~ 47 cm/s and c' ~ 12 cm/s, respectively, which 
are consistent with the values observed in human ventri- 
cles [13,16]. 

The reentrant wave activates the point in the ring 
(proximal to the zone of slow conduction) chosen to be 
the origin (x = 0), where the stimulus enters the ring, at 
time t — To. At time t = Ti, an activation wave is initi- 
ated through stimulation at x = 0. If this first stimulus 
is unable to terminate the reentry, a second stimulus is 
applied at t = T2, again at x = 0. Note that, the first 
stimulus is always able to terminate the reentry if it is 
applied when the region on one side of it is still refractory 
- leading to unidirectional propagation. This is identi- 
cal to the mechanism studied previously for terminating 
reentry by pacing within a ID ring [7]. However, in this 
study we are interested in the effect of pacing from a site 
away from the reentry circuit. In that case, it is generally 
not possible for the first stimulus to arrive at the reen- 
try circuit exactly at the refractory end of the reentrant 
wave (as discussed above). Therefore, we have used val- 
ues of Ti for which the first stimulus can give rise to both 
the anterograde, as well as the retrograde branches, and 
only consider reentry termination through block of the 
anterograde branch of the stimulated wave in the zone of 
slow conduction. 

Fig. 2 (top) shows an instance of successful ter- 
mination of the reentrant wave where the anterograde 
branch of the stimulus applied at T\ — 1146.22 ms 
(Ti - T = 428.00 ms) is blocked at the boundary of the 
zone of slow conduction (at time t ~ 1253 ms, x = 50 
mm). Fig. 2 (bottom) shows a magnified view of the 
region at which conduction block occurs. The region 
immediately within the zone of slow conduction shows 
depolarization but not an action potential (i.e., the ex- 
citation is subthreshold), while the region just preced- 
ing the inhomogeneity has successively decreasing action 
potential durations. The peak membrane potential at- 
tained during this subthreshold depolarization sharply 
decreases with increasing depth into the inhomogeneity 
so that, at a distance of 0.5 mm from the boundary in- 
side the zone of slow conduction, no appreciable change 
is observed in the membrane potential V. For the same 
simulation parameters as Fig. 2, a single stimulus applied 
at T\—Tq > 428.63 ms is not blocked at the inhomogene- 
ity, and a second stimulus needs to be applied at time T 2 
to terminate reentry. 

Different values of coupling interval (T\ — To) and pac- 
ing interval (T 2 — Ti) were used to find which parameters 
led to block of the anterograde wave. Fig. 3 is a param- 
eter space diagram which shows the different parameter 



regimes where termination was achieved. 

The conduction block of the anterograde branch of the 
stimulated wave is due to a dynamical effect linked to a 
local increase of the refractory period at the boundary 
of the zone of slow conduction in the reentrant circuit. 
To obtain an idea about the variation of refractory period 
around the inhomogeneity, we measure the closely related 
quantity, action potential duration (APD). Fig. 4 shows 
the variation of the APD as the reentrant wave propa- 
gates across the ring. The APD of an activation wave 
was approximated as the time-interval between succes- 
sive crossing of —60 mV by the transmembrane potential 
V [17]. As the wave crosses the boundary into the region 
of slow conduction, the APD increases sharply. Within 
this region, however, the APD again decreases sharply. 
Away from the boundaries, in the interior of the region 
of slow conduction, the APD remains constant until the 
wave crosses over into the region of faster conduction 
again, with a corresponding sharp decrease followed by a 
sharp increase of the APD, around the boundary. This 
result is consistent with the observation by Keener [18] 
that non-uniform diffusion has a large effect on the re- 
fractory period in discrete systems. 

However, the APD lengthening alone cannot explain 
the conduction block. Fig. 4 shows that the maximum 
APD, at the boundary of the inhomogeneity, is ~ 345 ms, 
yet conduction failed for a coupling interval of 428.00 ms. 
The blocking of such a wave (i.e., one that is initiated 
at a coupling interval that is significantly longer than 
the previous APD) implies that APD actually underes- 
timates the effective refractory period. While it is be- 
lieved that such "postrepolarization refractoriness" does 
not occur in homogeneous tissue, it has been observed 
experimentally for discontinuous propagation of excita- 
tion [19,20]. Postrepolarization refractoriness is caused 
by residual repolarizing ionic current that persists long 
after the transmembrane potential has returned almost 
to its quiescent value. As seen in Fig. 5, this has oc- 
curred in our system - at the time when the electrotonic 
current of the stimulated wave begins entering the inho- 
mogeneity, the ionic current has not fully recovered from 
the previous wave, and is actually working against the 
electrotonic current (7 io „ > 0). If the stimulated wave is 
early enough, this ionic current mediated postrepolariza- 
tion refractoriness is enough to prevent the electrotonic 
current from depolarizing the zone of slow conduction. 

This prolongation of the refractory period at the bor- 
der of the inhomogeneity is the underlying cause of the 
conduction block leading to successful termination of 
reentry, as can be seen in Fig. 6. Here, we compare 
the behavior at the inhomogeneity boundary when the 
wave is blocked (CI = 428.00 ms) with the case when the 
wave does propagate through the inhomogeneity (CI = 
428.64 ms). As already mentioned, in the case of block, 
the depolarization of the region immediately inside the 
zone of slow conduction is not sufficient to generate an ac- 
tion potential. Instead, the electrotonic current increases 
the transmembrane potential V to <~ — 55 mV and then 
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slowly decays back to the resting state. On the other 
hand, if the stimulation is applied only 0.64 ms later, V 
(as in the earlier case) rises to ~ —55 mV and then after 
a short delay (~ 8 ms) exhibits an action potential. This 
illustrates that in the earlier case the region had not yet 
recovered fully when the stimulation had arrived. 

Fig. 7 shows in detail how the difference in diffusion 
constant in the inhomogcneity leads to conduction block 
[21]. In the region with normal diffusion constant [Fig. 
7 (a)], the initial stimulation for a region to exhibit an 
action potential is provided by current arriving from a 
neighboring excited region. This spatial "electrotonic" 
current I spat — C m D\7 2 V is communicated through gap 
junctions which connect neighboring cardiac cells. As 
soon as the "upstream" neighboring region is excited with 
a corresponding increase in potential V, the difference in 
membrane potentials causes an electrotonic current to 
flow into the region under consideration (positive deflec- 
tion in Ispat)- This causes the local potential to rise and 
subsequently, an outward current to the "downstream" 
neighboring non-excited region is initiated (negative de- 
flection in I sp at)- Thus, the initial net current flow into 
the region is quickly balanced and is followed, for a short 
time, by a net current flow out of the region, until the 
net current Ispat becomes zero as the membrane poten- 
tial V reaches its peak value. The change in V due to 
Ispat initiates changes in the local ionic current Ii on . This 
initially shows a small positive hump, followed by a large 
negative dip (the inward excitatory rush of ionic cur- 
rent), which then again sharply rises to a small positive 
value (repolarizing current) and gradually goes to zero. 
The net effect of the changes in I spat and I ion on the 
local membrane potential V is reflected in the curve for 
dV/dt = DV 2 V - {h on /C m ). This shows a rapid in- 
crease to a large positive value and then a decrease to a 
very small negative value as the membrane potential V 
reaches its peak value followed by a slight decrease to the 
plateau phase value of the action potential. 

At the boundary of the inhomogcneity, however, the 
current I spa t is reduced drastically because of the low 
value of the diffusion constant D' at the inhomogeneity 
[Figs. 6(b) and 7(b)]. Physically, this means that the net 
current flow into the cell immediately inside the region of 
slow conduction is much lower than in the normal tissue. 
This results in a slower than normal depolarization of 
the region within the inhomogcneity. Because the region 
is not yet fully recovered, this reduced electrotonic cur- 
rent is not sufficient to depolarize the membrane beyond 
the excitation threshold. Therefore, the depolarization is 
not accompanied by the changes in ionic current dynam- 
ics needed to generate an action potential in the tissue 
[Fig. 6(c)]. The resulting dV/dt curve therefore shows 
no significant positive peak. Hence, the initial rise in V 
is then followed by a decline back to the resting potential 
value as is reflected in the curve for dV/dt in Fig. 6(d). 

When the stimulation is given after a longer coupling 
interval (T\ — T ) , i.e., after the cells in the inhomo- 
geneous region have recovered more fully, a different se- 



quence of events occurs. As in the case shown in Fig. 
7(b), the depolarization of the region immediately inside 
the zone of slow conduction is extremely slow because of 
the reduced electrotonic current I spa t- However, because 
the tissue has had more time to recover, the electrotonic 
current depolarizes the membrane potential beyond the 
threshold and the ionic current mechanism responsible 
for generating the action potential is initiated (as illus- 
trated by the broken curves in Fig 6). As a result, the 
excitation is not blocked but propagates through the in- 
homogeneity, although with a slower conduction velocity 
than normal because of the longer time required for the 
cells to be depolarized beyond the excitation threshold. 

To ensure that these results are not model dependent, 
especially on the details of ionic currents, we also looked 
at a modified Fitzhugh-Nagumo type excitable media 
model of ventricular activation proposed by Panfilov [22]. 
The details about the simulation of the Panfilov model 
are identical to those given in Refs. [23,24]. This two- 
variable model lacks any description of ionic currents 
and does not exhibit either the restitution or dispersion 
property of cardiac tissue. As expected, such differences 
resulted in quantitative changes in termination require- 
ments (e.g., the parameter regions at which termination 
occurs for the Panfilov model are different from those 
shown in Fig. 3 for the Luo-Rudy model). Nevertheless, 
despite quantitative differences, these fundamentally dif- 
ferent models shared the requirement of an inhomogene- 
ity for termination, thereby supporting the model inde- 
pendence of our findings. 

There are some limitations of our study. The most 
significant one is the use of a ID model. However, our 
preliminary studies on pacing in a 2D excitable media 
model of anatomical reentry [24] show similar results. 
We have also assumed the heart to be a monodomain 
rather than a bidomain (which has separate equations 
for intracellular and extracellular space). We believe this 
simplification to be justified for the low antitachycardia 
pacing stimulus amplitude. We have used a higher degree 
of cellular uncoupling to simulate ischemic tissue where 
slow conduction occurs. However, there are other ways 
of simulating ischemia [25], e.g., by increasing the exter- 
nal K + ion concentration [26]. But this is transient and 
does not provide a chronic substrate for inducing irregu- 
lar cardiac activity. Also, there are other types of inho- 
mogeneity in addition to ischemic tissue. For example, 
existence of a region having longer refractory period will 
lead to the development of patches of refractory zones 
in the wake of the reentrant wave. If the anterograde 
branch of the stimulus arrives at such a zone before it 
has fully recovered, it will be blocked [1]. 

Despite these limitations, the results presented here of- 
fer insight into pacing termination of anatomical reentry 
in the ventricle. We have developed general (i.e., model 
independent) mathematical arguments, supported by 
simulations, that circuit inhomogeneities are required for 
successful termination of anatomical-reentry VT when 
stimulation occurs off the circuit (as is typical in real- 
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ity). Thus, considering the critical role of such inhomo- 
geneities may lead to more effective pacing algorithms. 
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FIGURES 




(c) (d) 

FIG. 1. Schematic diagram of reentry in a ID ring illustrat- 
ing the necessity of a region of inhomogeneity for successful 
termination of reentry by pacing. At t — the reentrant wave 
reaches the entrance sidebranch (a). At t — z/2c the reentrant 
wave propagating through the sidebranch and the first stim- 
ulus mutually annihilate each other (b). At t = r + (z/c) the 
second stimulus reaches the reentry circuit by which time the 
refractory tail is a distance z away from the sidebranch (c). 
The presence of a region of inhomogeneity in the ring makes 
it possible that the anterograde branch of the second stimulus 
will encounter a refractory region behind the reentrant wave 
(d). 
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FIG. 2. (top) Plot of the membrane potential, V , show- 
ing spatiotemporal propagation of a reentrant wave in a 
Luo-Rudy ring of length 250 mm, successfully terminated by 
pacing with a single stimulus. The zone of slow conduction is 
between x = 50 mm and x = 75 mm. In this region the dif- 
fusion constant changes from D = 0.556 to D' — 0.061 cm 2 /s 
with an infinite gradient at the boundaries. The reentrant 
wave activates the site at x = mm at To = 718.22 ms. The 
stimulus is applied at x = mm at T\ — 1146.22 ms (coupling 
interval = 428.00 ms). (bottom) Magnification of the above 
plot for the region between x = 43 mm and x — 51 mm and 
the time interval 1242.50 < t < 1277.50 ms. Note that, inside 
the inhomogeneity, no action potential is generated. 
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FIG. 3. Parameter space diagram of coupling interval (CI) 
and pacing interval (PI) at which termination occurs in the 
ID Luo-Rudy ring of length 250 mm with a zone of slow con- 
duction between x = 50 mm and 75 mm (D = 0.556 cm 2 /s, 
D' = 0.061 cm 2 /s). The VT period around the ring is 692.46 
ms. The line connecting the circles represents the critical 
pacing interval value below which the second stimulus gets 
blocked by the refractory tail of the first stimulus. The re- 
gion between the circles and squares represents the regime 
in which the second stimulus is blocked in the anterograde 
direction in the zone of slow conduction (leading to success- 
ful termination) for a boundary with infinite gradient. For 
longer PI (the region above the squares), the second stimu- 
lus propagates through the inhomogeneity and only resets the 
reentrant wave. Note that the APD of the reentrant wave at 
x = is approximately 340.6 ms. For CI < 355.29 ms, the 
first stimulus is blocked at x = 0, while for CI > 428.63 ms, 
the anterograde branch of the first stimulus is not blocked 
at the zone of slow conduction (leading to resetting of the 
reentrant wave). For intermediate values of CI (i.e., between 
355.29 — 428.63 ms), the first stimulus is able to successfully 
terminate reentry. 
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FIG. 4. Variation of the action potential duration (APD) 
for the reentrant wave as it proceeds from x — to x = 125 
mm in the ring. The arrow indicates the direction of wave 
propagation. The dashed lines between x = 50 mm and 75 
mm enclose the region where the diffusion constant changes 
from D — 0.556 to D' — 0.061 cm 2 /s with an infinite gradient 
at the boundaries. 
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FIG. 5. The ionic current Ii 0n (solid curve) and the spa- 
tial electrotonic current I spa t = C m DX7 2 V (broken curve) at 
x — 50 mm (the boundary of the inhomogeneity). Note that, 
initially, the latter is two orders of magnitude smaller than 
the former. At t ~ 1248.5 ms, I spa t begins to increase due 
to current inflow from the upstream excited neighboring re- 
gion. The high value of I ion throughout this time indicates 
that the boundary of the inhomogeneity has not yet recovered 
from the passage of the previous wave by the time the next 
stimulated wave arrives. 
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FIG. 6. The (a) transmembrane potential V, (b) the spatial 
(electrotonic) current I spa t = C m DV 2 V , (c) the ionic current 
lion, and (d) the rate of change of transmembrane potential 
dV/dt = DV 2 V - {lion /C m ), at x = 50 mm for coupling 
intervals CI = 428.00 ms (solid curves) and 428.64 ms (broken 
curves) [D — 0.556 cm 2 /s, D' — 0.061 cm 2 /s]. In the former 
case, the initial depolarization is insufficient to generate an 
action potential and the excitation wavefront is blocked. In 
the latter case, the region has recovered sufficiently so that 
the stimulation is able to generate an action potential and the 
wave propagates through the inhomogeneity. 
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FIG. 7. The spatial variation of transmembrane poten- 
tial V, the spatial current I spa t = C m D\7 2 V, the ionic cur- 
rent hon, and the rate of change of transmembrane potential 
dV/dt = DX? 2 V — (lion I Cm), in (a) a homogeneous region 
of tissue (at t = 1240.05 ms) and (b) at the border of slow 
conduction (at t = 1254.00 ms). The zone of slow conduction 
is between x = 50 mm and x = 75 mm (D = 0.556 cm 2 /s, 
D' — 0.061 cm 2 /s) and the coupling interval is 428.00 ms. In 
(a) the excitation wavefront propagates normally, while in (b) 
the wavefront is blocked at the boundary of the inhomogene- 
ity. The arrows indicate the direction of propagation of the 
wavefront. 
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